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^ ! ABSTRACT 

■ The Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) provides an 
• • ' unprecedented opportunity to search for blazars at sub-mm wavelengths. We cross- 

, matched the FIRST radio source catalogue with the 11655 sources brighter than 35 

k> ■ mJy at 500 ^m in the ~ 135 square degrees of the sky covered by the H-ATLAS 

\ equatorial fields at 9 h and 15 h, plus half of the field at 12 h. We found that 379 of the 

C$ • H-ATLAS sources have a FIRST counterpart within 10 arcsec, including 8 catalogued 

blazars (plus one known blazar that was found at the edge of one the H-ATLAS 
maps). To search for additional blazar candidates we have devised new diagnostic 
diagrams and found that known blazars occupy a region of the log^soo/jm/Saso/mi) 
vs. log(55oopm/5'i.4GHz) plane separated from that of the other sub-mm sources with 
radio counterparts. Using this diagnostic we have selected 12 further candidates that 
turn out to be scattered in the (r — z) vs. (u — r) plane or in the WISE colour-colour 
diagram proposed by Massaro et al. (2012), where known blazars are concentrated in 
well defined strips. This suggests that the majority of them either are not blazars or 
have spectral energy distributions contaminated by their host galaxies. A significant 
fraction of true blazars are found to be hosted by star- forming galaxies. This finding, 
supported by an analysis of blazars detected in Planck 545 and 857 GHz bands, is 
at odds with the notion that blazar hosts are passive ellipticals and indicates that 
the sub-mm selection is providing a novel prospect on blazar properties. Based on 
an inspection of the available photometric data, including the WISE all-sky survey, 
the unpublished VIKING survey and new radio observations, we tentatively estimate 
that there are 11 blazars with synchrotron flux density S^oo^m > 35mJy over the 
considered area. This result already allows us to constrain blazar evolution models. 
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1 INTRODUCTION 

Blazars are a subclass of active galactic nuclei (AGNs) char- 
acterized by non-thermal continuum emission from radio to 
7-rays. They have extreme properties: they are strongly vari- 
able across the full electromagnetic spectrum, reach very 
high observed luminosities and are frequently strongly po- 
larized. Their spectral energy distribution (SED) is charac- 
terized by two broad peaks in uL v . The first peak occurs at a 
frequency Vp varyin g from ~ 10 12 Hz to ~ 10 19 Hz (~ 300/xm 
to ~ 3 x 10" s /im) l|Nieppola, Tornikoski fc Valtaoial [200^ 1 
and is attributed to Doppler boosted synchrotron emission 
from a highly relativistic jet pointing towards the observer. 
The continuum spectrum is "flat" (spectral index a > —0.5, 
S v oc u a ) from radio to sub-mm or even shorter wavelengths. 
The second peak, attributed to inverse Compton scattering, 
occurs at 7-ray energies. 

However, the non-thermal emission not always out- 
shine other components, like emission from the host galaxy, 
from the accretion disk and from the circum-nuclear torus, 
to the point of making such components undetectable. 
Blazar hosts are generally found to be early-type galaxies 
(iKotilainen. Falomo. fc Scarpalll998l;IO'Dowd fc Urry||2005l : 



iKotilainen et all 20071 ; Leon-Tavares et al. 2011 ) and may 
therefore contaminate the non-thermal spectrum primarily 
at optical/near-IR wavelengths. However, we should also 
take into account the possibility that at least some of 
them are endowed with star formation activity associated 
with thermal dust emission at far-IR/sub-mm wavelengths. 
The direct thermal emission from the accretion disk may 
show up as a relatively narrow 'bump' in the optical/UV 
iDermer fc Schlickeiserlll993l ; iGhisellini fc Tavecchid 120091 ') 
while the thermal emission from the dusty circum-nuclea r 
torus may be detectable in the mid-IR (|Perlman et al.f 2008). 

The optical spectra of blazars show a striking di- 
chotomy, leading to the recognition of two sub-classes: BL 
Lacertae objects (BL Lacs) with an almost featureless spec- 
trum, and Flat-Spectrum Ra dio Quasars (FSRQ s) with 
strong, broad emission lines. iPadovani fc Giommil (l995) 
classified the BL Lacs on the basis of the synchrotron peak 
frequency, t/p, into low, intermediate and hig h frequency 
peake d. This classification was extended by lAbdo et al.l 
|20ld ) to all blazars: low synchrotron peak (LSP: Vp < 
10 14 Hz), intermediate synchrotron peak (ISP: 10 14 < Vp < 
10 15 Hz) and high synchrotron peak (HSP: v% > 10 15 Hz) 
objects. The distribution of Vp is bimodal, with only a mi- 
nor fraction of objects peaking at intermediate frequencies. 
The bimodality however may be, at least partly, due to the 
blazar selection which, so far, was mostly done either in the 
radio or in the X-ray band. The former selection favours 
LSPs, the latter HSPs. 

The Herschel selection, at intermediate frequencies, 
should help us to pick up ISPs, thus providing constraints on 
the abundance of these objects and indications on whether 
or not there is a continuity in the blazar population, from 
LSPs to HSPs. Moreover, a full frequency coverage is crucial 
to disentangle the different, non-thermal and thermal, con- 
tributions to the blazar SEDs, as briefly summarized above. 
The Hersch el Astrophysical Terahertz Large Area Survey 
fH- ATLAS; lEales et al. 20101 ). the largest ar ea survey car- 



ried out by the Herschel space observatory JPilbratt et al 
2010) covering ~ 550 deg 2 with PACS (|Poglitsch et al 



l201ll ) and SPIRE (|Griffin et alj|2010h instruments between 
100 and 500 /im, provides an unprecedented opportunity 
to obtain flux density limited blazar samples at sub-mm 
wavelengths. As for the non-thermal emission, the Her- 
schel data cover a particularly interesting frequency range 
close to the synchr otron peak of the most luminous LSPs 
(Fossati et al. 1998),where information is currently scanty 
(jpadovani et al. IbOOfJ V Since Vp carries information on the 
Lorentz factor 7 P of emitting electrons, on the Doppler fac- 
tor 5 and on the magnetic field strength B (y v oc jpSB), 
its value provides crucially important information on key 
physical parameters. Moreover, the sub-mm selection is es- 
pecially well suited to look for signatures of thermal dust 
emission from the host galaxy, thus providing a new perspec- 
tive on blazar hosts and allowing us to test the early-type 
host paradigm. 

Blazars are however a tiny fraction of H-ATLAS 
sources: their estimated surface density at the H-ATLAS 5 a 
detection limits is ~ 0.1 percent or less of the surface density 
of dusty galaxies. Singling them out is not easy and appro- 
priate diagnostic tools must be devised. An investigation 
of the H-ATLAS Science De monstration Phase (SDP) field 
|Gonzalez-Nuevo et alj201Ch uncovered two blazars brighter 
than th e average 5a dete ction limit at 500 /J,m, Ssoo/am — 
44mJy (|Rigbv et al.ll201 if ), over an area of 14.4 deg 2 . In this 
paper we extend the analysis to an area almost an order of 
magnitude larger. This allows us to build the first statisti- 
cally meaningful sample of sub-mm selected blazars. 

The layout of this paper is the following. In Section[2] 
we describe the selection of candidate blazars, present the 
new diagnostic diagrams we have devised, and report on 
the follow-up observations aimed at assessing the nature of 
the selected candidates. In Section[3]we briefly discuss their 
SEDs, built by combining the H-ATLAS photometry with 
literature and unpublished VIKING data. In Section[4] we 
discuss constraints on sub-mm blazar counts and implica- 
tions for evolutionary models. Finally, Section[S] summarizes 
our main results. 

We adopt a flat cosmology with J1a = 0.734, Q m = 
0.266, and h = .71 as derived from WMAP 7-yr data 
jLarson et al.ll201ll l. 



2 BLAZAR CANDIDATES 

In this study we use the catalogue of SPIRE sources de- 
tected in the two H-ATLAS equatorial fields centred at 
about 9h and 15 h, plus half of that at 12 h, with areas 
of 53.25 deg 2 , 53.93 deg 2 , and 27.37 deg 2 , respectively, for a 
total of 134.55 deg 2 . In the Phase 1 data release [Hoyos et 
al. (in preparation), Valiante et al. (in preparation)] data 
for the full 12 h field will be included. The catalogue con- 
tains sources that have been detected in the noise-weighted 
PSF-filtered 250 fim map app lying the MADX algorithm 
|Maddox et al.. in preparation! - ) . Then, at the position of 
the sources, flux densities are estimated at each of the 
SPIRE bands. Finally, only sources with a signal-to-noise 
^ 5 at any of the SPIRE bands are included in the cat- 
alogue. T he H-ATLAS maps and data reduct i on ar e dis- 
cussed in Pascale et "all (|201lh and llbar et al.l d2010h. th e 
source catalogue creation is described in Rig bv et al.l 
The SPIRE beams have FWHM of 18.1, 24.8 and 35.2 arc- 
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Figure 1. Diagnostic colour-colour diagrams: Ssoonm/Si^GHz vs - S^oo^m / SssO/j-m (upper panel) and (r — z) vs. (it — r) (lower panel). In 
both panels the grey dots denote the H-ATLAS sources with Ssoo^im ^ 35 mjy and FIRST counterparts (only a subset of these sources 
have SDSS photometry and thus appear also in the lower panel), except for the 9 catalogued blazars that are identified by the larger 
black + signs. The grey + signs enclosed in red squares denote the new blazar candidates. Both catalogued and candidate H-ATLAS 
blazars are labelled with their ID numbers (see Tables [3] and |4j . The Planck ERCSC blazars detected both at 545 and at 857 GHz 
(see text) are represented by the black X signs. In the upper panel the thick dashed line shows the colour-colour track of a prototype 
star-forming galaxy, Arp 220. The errorbars in iog(5 , 500/jm/>Si.4GHz) ^ or ca t a l°g UC( i an d candidate blazars are ^ 0.1 and have not 
been plotted for the sake of clarity. In the lower panel, the light blue X signs correspond to BZCAT blazars with SDSS photometry. In this 
panel we also show, for comparison, the colour-colour tracks, as a function of redshift, of 4 passive elliptical templates of different ages 
(red lines) and of 3 sta r-forming galaxies ( blue lines). The SO, M82, QSO and Arp220 tracks were computed using the SEDs tabulated 
in the SWIRE library jPoIletta et aljfeoOTl l available at |http:/ /www.iasf-milano. inaf.it/~polletta/t emplat es/swire_templat es . html| The 
elliptical galaxy tracks were computed using GALSYTH, the web based interface for the GRASIL model (Silva ct al. 1998) available at 
http://galsynth.oand.inaf.it/galsynth/. The colours corresponding to some redshift values are indicated. 
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Figure 2. Positional offset between the 379 H-ATLAS sources 
with SsoOMm ^ 35 mjy and their FIRST counterpart within 10 
arcsec. 

sec for 250, 350 and 5 00 /im respectively (|Swinvard et al.l 
|2010| ; iRigbv et al1l201ll ). The catalogue we have used con- 
tains 66535 unique objects (26369 in the 9h field, 12686 
in the 12h field, and 27480 in the 15h field). Since blazars 
have red sub-mm colours they are more easily singled out 
at the longer wavelengths. We have therefore confined our- 
selves to sources with Ssoo/jm ^ 35mJy, corresponding to 
4cr detections at this wavelength. There are 11655 sources 
satisfying this criterion. Blazars are powerful radio sources. 
Thus the first step towards spotting them is to require that 
they are detected by radio surveys covering our fields. A 
cross correlation of these sources with the FIRST (Faint 
Images of the Radio Sky at Tw enty-cm) survey catalogue 
|Becker. White, fc Helfandlll995h with a 10 arcsec search ra- 
dius yields 379 matches, out of the ~ 11000 FIRST sources 
present in the ~ 135 deg 2 H-ATLAS equatorial fields. 

The sample made of the 379 sources selected as de- 
scribed above (Ssoopm ^ 35 mjy and FIRST counterpart 
within 10 arcsec) will be referred to as the parent sample. 
The ty pical rms positiona l error for our sources is 2.5 arcsec 
or less (|Rigbv et al.| [20H r1. and the FIRST positions are ac- 
curate to better than 0.5 arcsec for Si.4GHz > 3mJy, so that 
the rms offset between the radio and the sub-mm position 
is a ^ 2.55 arcsec (see Fig. [2]). In the case of a Gaussian 
distribution of positional errors the probability that a true 
counterpart has an apparent positional offset ^ A is 

p(> A) = exp[-0.5(A/(r) 2 ]. (1) 

For A = 10 arcsec, p(> A) ^ 4.6 x 10~ 4 and the number of 
true FIRST counterparts that we may have missed with our 
11655 trials is ~ 5, i.e. ~ 1.5 percent of the parent sample. 
The surface density of FIRST sources is wfirst — 

The quoted positional error refers to > 5 a 250 fim sources as 
most (355) of the 379 sources in the parent sample are. 



Blazar population 



Obscured AGNs 
•16 




i ... i .... i .... ' .... i .... i 
1 2 3 4 5 

[4.6]-[12] 



Figure 3. Diagnostic diagram using WISE colours. BZCAT ob- 
jects are shown in cyan and our candidate and confirmed blazars 
are labelled with their ID numbers (see Tables [3] and 2). The 
other da ta points and the boun daries of the 'WISE blazar strip' 
are from lMassaro et al.l j2011ah . BZB, BZQ, and BZU stand for 
BL Lac, fiat-spectrum radio quasar, and blazar of uncertain type, 
respectively. The contours tracing the density of sources in colour 
bins with annotations showing the location of differe nt classes of 
objec ts are from the WISE Explanatory Supplement iCutri et al.l 
l201ll) . 

90 deg~ 2 . Therefore the probability that a FIRST source lies 
by chance within A = 10 arcsec from a given Herschel source 
is 7rA 2 npiRST — 2 x 10~ 3 and the expected number of chance 
associations is ~ 25. However, as we will see in the following, 
from the distribution of Ssoopm/ SiagHz flux density ratios 
of known blazars we infer that genuine blazars must have 
S , 5Q0/im/5i.4GH« < 2.4 or, in our case, Si.4GHz > 14mJy. 
The surface density above this limit is ~ 13.5 deg -2 and the 
number of chance associations decreases to ~ 4. 

2.1 Diagnostic diagrams 

A characteristic property of blazars is their peculiar contin- 
uum spectrum. This property allows us to construct diag- 
nostic diagrams useful to select the tiny fraction of blazar 
candidates out of the overwhelming number of dusty galax- 
ies. The relatively flat non-thermal continuum of blazars 
from radio to sub-mm wavelengths implies that they have 
red sub-mm colours and sub-mm to radio flux density ratios 
substantially lower than those of dusty galaxies. Thus we 
expect that blazars occupy in, e.g., the log(S , 5oo / jm/5 , 350Mm) 
vs log(S , 500Mm/5'i.4GHz) plane a region separated from that 
populated by sub-mm galaxies dominated by thermal dust 
emission. 

This expectation is borne out by the upper panel of 
Fig. [T] To locate the "blazar region" in this diagram we 
have looked for known blazars among the 379 sources in the 
parent sample by cross-matching them wit h the most re- 
cent v ersion of the blazar catalogue BZCAT (|Massaro et al.l 
l2011bl ). again with a search radius of 10 arcsec. We got 
8 matches (see Table |3J. Since the BZCAT positions are 
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mostly taken from radio catalogues and are generally bet- 
ter than 1 arcsec, after eq. ^ the expected number of real 
matches that may have been missed with our 379 trials is 
^ 0.4. The probability that a BZCAT source lies by chance 
within 10 arcsec from a given Herschel source is difficult to 
assess accurately because of the inhomogeneous distribution 
of BZCAT sources. However the BZCAT catalogue contains 
~ 3100 sources over the whole sky so that their mean surface 
density is wbzcat < 0.1 deg" 2 . If they were uniformly dis- 
tributed the expected number of chance associations would 
be < 10 -3 . There is therefore a wide enough margin to con- 
fidently assume that all the associations found are real even 
though the distribution is not uniform. In particular, of the 
8 matches found with a search radius of 10 arcsec, the two 
with the largest angular distance to a FIRST source (ob- 
jects #1 and 2 of Table [3l with a separation of 8.67 and 
9.72 arcsec, respectively) have Si.4GHz > 300 mjy. Since the 
surface density of FIRST sources brighter than 300 mjy is 
~ 0.46 deg -2 , the probability of a chance association within 
10 arcsec is ~ 1.1 x 10 -5 and the expected number of chance 
associations within 10 arcsec is 0.13. On the other hand, af- 
ter eq. Q the probability that a true counterpart has a 
positional offset of 9.72 arcsec or of 8.67 arcsec is ~ 7 x 10 -4 
or ~ 3 x 10 -3 , respectively. Thus the blazar classification of 
these 2 H-ATLAS sources must be taken with caution. 

In the BZCAT catalogue there are an additional 15 
blazars located in our fields. For only two of them could we 
detect a significant signal in the H-ATLAS noise-weighted 
beam convolved maps. One (BZBJ0849+0206) has a 500 
flux density of 31.6 mjy (S/N ~ 3.6) and is weaker at 
250 and 350 /im (S/N ~ 2.1 and 2.2, respectively). Since 
this object is fainter than the adopted flux-density limit 
(Ssooum = 35mJy) we have not considered it further. The 
other object (BZQJ1404-0013) has Ssoo^m = 45 mjy, i.e. is 
a 5cr detection at this wavelength, but is below the 5a de- 
tection limit at 250 and 350 /tm (^so/jm — 18 ± 7mJy and 
S350/jm — 32 ± 8mJy). It was probably missed because of 
its low 250 fim flux density (cand idate sources were identi- 
fied as > 2.5cr peaks at 250 /im; Ri gbv et al.ll201lh ; more- 
over, it lies close to the border of the field. We have added 
this source to our blazar sam ple. We have also searched the 
2-vr FermirLAT catalogu<EI (|The Fermi-LAT Collaboration! 
l201ll ) for counterparts of our sources that might be blazars 
not listed in the BZCAT. None was found: the only Fermi- 
LAT source located within our fields and not listed in the 
BZCAT is identified as a pulsar and is about 3 arcmin away 
from the nearest source in the parent sample. 

To increase the sample of known blazars with measured 
sub-mm colours we have cross-correlated with the BZCAT 
the full-sk y Planck Early Release Compa ct Source Catalogue 
(ERCSC; IPlanck Collaboration! boilal ). which consists of 
sources with flux density measurements at both 545 and 857 
GHz. We found 14 matches The additional data available in 
the NASA Extragalactic Databas<0 (NED) show that one of 
these sources, BZUJ1325-4301, is the well known extended 
radio galaxy Centaurus A, while the blazar classification 
of BZBJ1 136+1601 is dubious and BZUJ2209-4710 (alias 
NGC 7213) has sub-mm colours indicative of a strong con- 



2 fcrmi .gsfc. nasa.gov /ssc /data/ access /lat / 2yr .catalog/ 
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tamination by thermal dust emission from the host galaxy. 
Moreover, the ERCSC flux densities of BZQJ1559+0304 and 
BZQJ1719+0817 may be affected by source blending. The 
flux density ratios for these sources are therefore not rep- 
resentative of the blazar non-thermal continuum. We are 
then left with 9 blazars (see Table [5}, whose colours are 
also plotted in the upper panel of Fig. [1] According to the 
technical specifications of their respective instruments, the 
Planck/RFI 857 GHz and the Herschel/ SPIKE 350 urn chan- 
nels have almost exactly the same central wavelength and 
roughly the same bandwidth. The comparison between the 
Planck 545 GHz and the SPIRE 500 flux densities re- 
quires a colour correction because the central frequency of 
the 500 /im channel is somewhat higher (~ 600 GHz). The 
extrapolation from 545 GHz to 600 GHz was made using, for 
each Planck blazar, its spectral index between 545 and 857 
GHz; since the spectral indices are quite flat, the corrections 
are small. 

The upper panel of Fig. [1] shows that the 9 known 
blazars in our H-ATLAS fields and the 9 ERCSC 
blazars lie in a region in the log(S , 5 oo,im/>S'i.4GHz) vs 
log(S , 5 ooMm/S , 35o / jm) plane fairly well separated from the re- 
gion populated by the majority of the other sources in our 
parent sample that have a Ssoo/jm/Si^GHz ^ 10, typical of 
star forming galaxies. To illustrate this further, we have in- 
cluded in this panel a thick dashed line that represents the 
colour track as a function of redshift of the prototype star- 
burst galaxy Arp 220. Objects with <S , 5oo Mm /S'i.4GHz < 10 
likely have radio emission dominated by an AGN compo- 
nent, that since we are using low-frequency (1.4 GHz) radio 
observations, is expected to be mostly steep-spectrum (i.e. 
not associated to blazars). 

The 5 , 5oo^m/Si.4GHz ratios of blazars are ^ 2.4 and their 
SsooMm/'S'asoMm colours are redder than average. This upper 
limit to S , 500fjm/S'i.4GHz ratio for blazars comes from the 
ERCSC objects, while the H-ATLAS ones, which may be 
thought to be more relevant here, have significantly lower 
ratios. But we must beware of a selection effect. Since the 
BZCAT blazars are mostly radio selected (with important 
additions from X-ray and 7-ray selections), relatively radio- 
faint objects, such as possible H-ATLAS blazars with high 
Ssooum/ SiAGHz ratios, are easily missed by previous blazar 
searches. In fact, as will be seen in the following, our candi- 
date H-ATLAS blazars have very limited photometric data. 
The situation is different for the much (generally more than 
10 times) brighter ERCSC blazars for which we likely have 
the full range of Ssooum/ Siaghz ratios. We have there- 
fore adopted as the limiting ratio, the maximum value for 
ERCSC blazars. Indeed the main strength of the sub-mm 
selection compared to the traditional ones is its unique effi- 
ciency in selecting objects with high Ssoo/jm/Si^GHz ratios. 
This means that even with the adopted limiting ratio we 
may miss some blazars. For this reason we plan follow-up 
observations also of objects with ratios above the chosen 
limit although we expect that such objects are unlikely to 
be blazars. 

The criterion SsooMm/Si^GHz ^ 2.4 yields 12 new candi- 
date blazars in our parent sample (see Table [4}. From them, 
7 have an angular distance to a FIRST source > 5 arcsec and 
a Si.4GHz > 23 mjy. Since the surface density of the FIRST 
sources at this flux density limit is ~ 8 deg -2 , the expected 
number of chance associations for angular distances between 
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5 and 10 arcsec is ~ 1.7 and some of these associations may 
indeed be spurious. Thus almost of all the objects in our 
parent sample are likely to be pure star-formin g galaxies or 
radio galax ies in star-forming hosts as studied in ljarvis et al.l 
l|2010l ) and lHardcastle et al.l l|2010h respectively. 

As a second diagnostic tool, we have built the (r — z) 
vs. [u — r) colour-colour diagram (lower panel of Fig. [1]) 
for kno wn blazars by cross -correlating the SDSS DRqj cat- 
alogue l|Aihara et al.ll201lh with the BZCAT, using a search 
radius of 0.5 arcsec. The surface density of DR8 sources is 
~ 3.22 x 10 4 deg -2 so that the probability of a chance as- 
sociation within 0.5 arcsec is ~ 2 x 10 -3 . Since there are 
~ 1800 BZCAT sources within the DR8 area, we expect 
~ 3.5 chance associations. We get 1600 matches. They are 
mostly concentrated along a relatively narrow strip in the 
(l — z) vs. (u — r) plane. Seven out of the 12 new candidates 
(# 11, 15, 16, 17, 19, 20, and 21) lie within or very close 
to this strip, although the Herschel/ SPIRE colours of # 16 
and 17 may be indicative of contamination by dust emission 
in the host galaxy. Objects # 12, 14, and 18 are at, or some- 
what above, the red end of the strip, suggesting that their 
optical colours may be moderately contaminated by the host 
galaxy. Hints of host galaxy contamination come also from 
the Herschel/ SPIRE colours of two of these objects (# 12 
and 14). As illustrated by the colour tracks in the lower panel 
of Fig. [1] very red colours are indicative of contamination 
by a passive early-type host galaxy. Indications that colours 
redder than (u—r) = 1.4 lik ely imply host galaxy contamina- 
tion w ere also reported bv lMassaro. Nesci. fc Piranomontg 
l|2012j ). The sub-mm to optical photometric data on one of 
the two outliers with blue (u — r) colours (object # 10) 
are strongly indicative of a star-forming galaxy. The other 
blue outlier (object # 13) is anomalously faint in the r-band 
(compared to neighbour bands); its (r — z) versus (u — r) 
colours may thus not be a good diagnostic. Some of the new 
Blazar candiates have optical colours consistent with QSOs. 

A third diagnostic tool, based on the [3.4] [4.6] [121 M m 
colour-colour diagram using WISE JWright et all l201Ch 
magnitudes, was proposed bv lMassaro et alT ( 2011al ). Mag- 
nitudes for the 3 WISE channels are availably for 7 H- 
ATLAS catalogued blazars (the exceptions are # 7 and 9) 
and for 10 H-ATLAS candidates (the exceptions are # 20 
and 24). However objects # 9 and 20 had magnitudes listed 
in the Preliminary Release of WISE data. The positions of 
H-ATLAS confirmed or candidate blazars in this diagram 
are identified, in Fig. [3] by their ID numbers (see Tables [3] 
and©. The catalogued H-ATLAS blazars for which WISE 
photometry is available plus the candidates # 10, 16, and 
20 lie in the upper part of 'blazar strip', i.e. in the region 
occupied by fiat-spectrum radio quasars (although # 16 is 
at the border). Four of our blazar candidates (# 12, 13, 
18, and 21) have WISE colours typical of starburst galax- 
ies, and two (# 17 and 19) are in the region populated by 
LINERs/obscured AGNfl 



ID RA (°) Dec. (°) S 1A GH* S 5 GHz °S 



10 


218.522 


-0.385 


59.92 


43.03 


5.25 


11 


220.332 


0.420 


57.32 


38.72 


3.19 


13 


222.376 


2.610 


50.09 


46.93 


4.30 


14 


174.092 


0.815 


151.74 


37.21 


3.93 


21 


222.406 


0.693 


51.93 


55.20 


9.09 



4 http://www.sdss3.org/dr8/ 

5 http:/ /irsa. ipac.caltech.edu/Missions/wisc. html 

6 http:/ /wise2. ipac.caltech.edu/docs/releasc/prelim/expsup/fig- 
ures / sec2_2fl6_annot.gif 



Table 1. Flux densities (mjy) at 5 GHz and their errors (erg) 
measured with the Medicina radio telescope for 5 blazar candi- 
dates. The FIRST 1.4 GHz flux densities (mjy) are also shown, 
for comparison. 



2.2 Follow-up observations in radio 

We have performed follow-up observations with the Medici- 
na 32m single-dish telescope at 5 GHz of 6 blazar candidates. 
Observations were carried out between August 1st -16th 
2011 i n the OTF sca n mode IjMangum. Emerson, fc Greisenl 
120071 ; iRighini 1 120081 ') applying the observing strategy, the 
data reduction procedur es and the software to ols developed 
for the SiMPlE project (jProcopio et al.ll201lf ). Twenty-five 
hours were allocated for 5 GHz observations with a single 
feed receiver and a 150 MHz-wide band. The 37 arcmin 
length of each scan (~ 5 x HPBW) was covered at 3 ar- 
cmin/s in 12.5 s. A sampling rate of 40 ms allowed us to 
obtain 60 samples/beam. The scans were centred on the 
FIRST positions. 

Given the FIRST flux densities we do not expect to see 
any of our targets in a single scan, which, in optimal ob- 
serving conditions, should hence appear as a linear baseline 
that corresponds to the off-source zero level of the signal. 
Along the scans, Gaussian noise shows up as amplitude fluc- 
tuations. However, cloudy weather, the presence of random, 
but unfortunately not rare, contributions by radio frequency 
interference (RFI) or digital noise heavily affected portions 
of the data. These disturbances give rise to bumpy baselines 
or spike- like features. All the scans that do not show a linear 
profile or show such features were removed before running 
the data reduction pipeline. 

We have use d the OTF Scan Ca li bration-Reduc tion 
pipeline (OSCaR, iProcopio etldl l201ll ; IRighini I l201ll ) to 
get the conversion factor from raw data counts measured 
on calibrators to their known flux densities, to determine 
the component of the flux density error due to calibration, 
and to calculate the source flux densities and their errors. 
3C295 and 3C286 were used as flux density calibrators. One 
out of the six targeted sources was not detected because of 
technical problems. The 5 GHz flux densities for the other 
5 sources were measured with signal to noise ratios greater 
than 6. The results are given in Table[T] Source # 14 turned 
out to have a steep 1.4-5 GHz spectral index and is optically 
identified as a galaxy. For the other 4 sources (# 10, 11, 13, 
and 21) our observations indicate a flat spectral index, con- 
sistent with them being blazars. We caution, however, that 
the low frequency spectral index is not necessarily a con- 
clusive indication for or against a blazar classification since 
the low frequency radio emission frequently comes from a 
different component. Additional flux density measurements 
up to mm wavelengths are necessary to assess the nature of 
candidate blazars. Such observations are being planned. 
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Figure 4. SEDs of the 9 catalogued blazars in the H- ATLAS equatorial fields (Table [3L The red circles, magenta diamonds, blue 
squares, red 6-pointed stars, blue 5-pointed stars and green asterisks are H-ATLAS, WISE, VIKING, GMRT, FIRST and Medicina data, 
respectively. The other data points are from the NASA/IPAC Extragalactic Database (NED). The downward-pointing triangles are the 
results of subtracting from the data the fitted second order polynomial (solid lines) rep resenting the synch rotron emission. In panels 1 
and 6 the dot-dashed lines show the dusty galaxy template for log(L^ ust / Lq) > 11.5 (Smith et al. 2 0121 ). put at the blazar redshifts 
(z = 1.123 and z = 1.173) respectively. This template is only shown for illustration purposes and no fit to the data has been attempted: 
since the synchrotron component is highly variable and the data are not simultaneous, fits would have little meaning. The dashed curves 
on the right side of panels 3, 4, 6, 7, and 9 are black-body spectra representing the possible emission from the accretion disk. Again, 
these curves are not meant to be fits of the data. 



3 SPECTRAL ENERGY DISTRIBUTIONS 

Four out of the 9 known blazars in the H-ATLAS 
fields have counterparts in the 2-yr Fermi-LAT catalogue 
l|The Fermi-LAT Collaboration! |201ll . see Table |3j. The 
blazar J090910. 1+012135 happened to be observed twice by 
Herschel, on 2009 Nov 22 and 2010 May 30, and showed a 
strong variability, thus providing a direct confirmation of the 
non-thermal, jet-related nature of its sub-mm emission. Its 
250, 350 and 500 /im flux densities increased from 160, 194 
and 266 mjy, respectively, at the first epoch, to 258, 338 and 
425 mjy at the second epoch. This source was also found to 
be variable also in 7-rays: its Fermi-LAT variability index is 
159.5. 

Figure|3]shows the spectral energy distributions (SEDs) 
from radio to UV frequencies of the 9 known blazars 
in the H-ATLAS equatorial fields. To the data from the 



NASA Extragalactic Database (NED) and to the H-ATLAS 
data we have added the near-IR photometry provided by 
the VISTA Kilo-degree INfrared Galaxy sur vey (VIKING; 
ISutherland et ail I2OI2I; iFleuren et all [20121 ) and by the 
WISE surve\l1 l|Wright et ail l201fj| ). In addition, 5 of our 
catalogued blazars and 9 of our blazar candidates were de- 
tected by the Giant Metrewave Radio Telescope (GMRT) 
survey of H-ATLAS equatorial fields at 325 MHz (Mauch et 
al., in preparation). The measured flux densities are given 
in Tableland shown in Figs. H and [5] The 0.325-1.4 GHz 
(FIRST) spectral indices of 4 catalogued blazars are flat or 
inverted (0^325 > —0.5, S u oc v a ) while all the 9 blazar 
candidates have a-J^s < —0.5. Three of the latter objects 
(# 10, 15, and 18) have 1.4 GHz NVSS flux densities larger 



http: / /irsa. ipac.caltech.edu/Missions/wise. html 
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Figure 5. SEDs of the 12 candidate blazars in the H- ATLAS equatorial fields (Table [4). The meaning of data points and of the solid 
and dashed lines is t he same as in Fig. [4] The dot-dashed lines are host galaxy templates: a Seyfert 2 g alaxy from the SW IRE library 
llPolletta et aLlfepOTl) for objects # 10, 16, and 17; a dusty galaxy template with £og(Z/ dust / Xq) > 11.5 ISmith et alj|2012h for objects 
# 11 and 18; the Arp 220 template from the SWIRE library for objects # 12, 14, and 19. As for Fig. [4] these templates are shown for 
illustration only: no fit of the data was attempted. 



than the FIRST ones by a factor > 2. Since FIRST has a 
much higher angular resolution than NVSS, this suggests 
that these sources are extended. As noted in Section [2~2l 
however, a steep low-frequency spectrum is an indication, 
but is not necessarily a proof, that the object is not a blazar 
since the low-frequency emission may come from a compo- 
nent different from the relativistic jet. In addition, all these 
data are far from simultaneous. Hence any fit should be dealt 



with great caution and can only be considered as purely in- 
dicative. 

As usual (e.g., iFossati et al.l 1 19981 ). we fit the syn- 
chrotron peak with a second order polynomial (i.e., a 
parabola). The estimated rest-frame synchrotron peak fre- 
quencies (in terms of vL v ) are in the range 12.8 ^ 

logOpcak.Hz) < 14.2, with a median log I'peak.median.Hz - 

13.4. These objects are thus mostly LSPs, but with an 
extension in the ISP region. In only one case (# 3, see 
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ID 


5325MHz [mjy] 


ID 


5325MHz [mjy] 


1 


500 ± 15 


12 


255 ± 13 


2 


157 ± 5 


15 


155 ± 5 


3 


986 ± 31 


16 


156 ± 6 


4 


4354 ± 124 


17 


72 ± 5 


8 


815 ± 26 


18 


210 ± 7 


10 


602 ± 21 


19 


84 ± 15 


11 


283 ± 9 


21 


127 ± 12 



Table 2. Flux densities at 325 MHz and their errors obtained 
with the Giant Metrewave Radio Telescope for 5 confirmed blazar 
and 9 new candidates. 



IGonzalez-Nuevo et all l|2010tl ) do we have enough high- 
energy data to allow a meaningful estimate of the shape of 
the inverse Compton peak (not shown in the corresponding 
panel of Fig. [4]). For this object, the Compton dominance, 
defined as the ratio between the inverse-Compton and syn- 
chrotron peak luminosities is v Pl icLvp,lc /vp,sLvp,s — 8.6. 
Although 3 other blazars have 7-ray measurements, they 
are missing X-ray data and therefore the rising part of the 
inverse Compton peak is essentially unconstrained. 

Blazar # 3 also shows indications of a UV excess 
that can be attributed to emission from the accretion dis k 
l|Dermer fc Schlickeiser]|l993l; [chisellini fc Tavecchiol | 2009h . 
as already pointed out by IGonzalez-Nuevo et al.l l|2010l) . 
Some hints of a UV bump can also be discerned on the 
SEDs of blazars # 4, 6, 7, and 9. Under the standard as- 
sumption that the accretion disk emission is a combination 
of black-bodies with temperatures de pending on the distance 
from t he central black-hole [eq. (1) of I Chisellini fc Tavecchiol 
(2009)] the black hole mass, Mbh, can be estimated as 



M= 



10 9 Mr. 



1.2 



T 



2 x 10 4 K 



10 45 erg/s 



(2) 



where r\ is the mass to light conversion efficiency, T is the 
maximum black-body temperature, and Ld is the total lu- 
minosity of the accretion disk. The lack of sufficient simul- 
taneous data hampers reliable estimates of the key quan- 
tities T and Ld, hence of Mbh and of the accretion rate 
Mbh = L d /(r]C 2 ). Tentative values are M b h ~ few x 10 9 M Q , 
Mbh ~ 0.2-1 Mq yr , and Eddington ratios Ld/Lpdd ~ 
0.2-1 x 10" 2 . 

The Herschel/ 'SPIRE colours of blazars # 1 and 6 seem 
to indicate a contamination by a star-forming host galaxy 
with Lpir ~ 10 12 Lq, i.e. a star-f o rmation rat e SFR ~ 100- 
lTOMoyr" 1 <|KennicuttJ 1 19981 ; lLapi et al.l l201ll ). The 
association of a blazar with a star-forming galaxy is at odds 
with the notion that blazar hosts are passive ellipticals, as 
enerally found for radio, X-ra y and 7- ray selected blazars 



Kotilainen. Falomo. fc Scarpa! Il998l; lO'Dowd fc Urr 



20051; iKotilainen et al.l 120071 ; iLeon-Tavares et all 1201 



Urryl 
l201ll ; 



Giommi et al. 20111 ). To check the indication that a signif- 



icant fraction of sub-mm selected blazars may indeed be 
associated with star-forming galaxies we have inspected the 
photometric data available in the NED for 14 catalogued 
blazars with counterparts in the Planck ERCSC at both 
545 and 857 GHz (see Section [2~I) . We have excluded 
BZUJ1325-4301, i.e. Centaurus A, because of the difficulty 
of removing the contribution of the extended emission 



and BZBJ1136+1601 because its classification as a blazar 
is dubious. For 5 (i.e. ~ 40 percent) of these objects 
(BZUJ0840+1312, BZQJ0921+6215, BZQJ1559+0304, 
BZQJ1719+0817, and BZUJ2219-4710) the sub-mm con- 
tinuum appears to be dominated by thermal dust emission. 
Somewhat surprisingly, 3 of the 5 Planck ERCSC objects 
(those labelled BZQ) and both H-ATLAS blazars with pos- 
sibly/likely star- forming hosts are classified in BZCAT as 
flat-spectrum radio quasars, although quasar-type blazars 
are expected to outshine the host galaxy. 

As for our blazar candidates, although the available, 
generally non-simultaneous, data are insufficient to draw 
firm conclusions, especially on account of the strong blazar 
variability, we may note that 4 of them (# 13, 15, 20, and 
21) seem to have SEDs approximately described by a second 
order polynomial (see Fig. [5j, the usual representation of 
the synchrotron peak, and thus consistent with being blazar 
SEDs. The SEDs of the other candidates seem to be strongly 
contaminated, or dominated, by emission from a dusty star- 
forming galaxy. Extensive follow-up observations are neces- 
sary to assess the nature of these sources. Multi-frequency 
radio observations, already planned, can establish whether 
or not our candidates have the "flat" high-frequency spec- 
trum typical of blazars. Variability and polarization are fur- 
ther indicators of the blazar nature. 



4 NUMBER COUNTS 

The H-ATLAS survey allows us to put the first meaningful 
constraints on 500 /im blazar counts, thus providing a test on 
evolutionary models and, in conse quence, on the underly ing 
physics. For example, the model bv lDe Zotti et al.l (|2005r ) as- 
sumes, for all blazars, a flat radio spectral index followed by 
a parabolic decline above a synchrotron peak frequency that 
increases with decreasing radio luminosity a ccording to the 
'blaz ar sequence' model (|Fossati et al.l ll998). On the other 
hand iTucci et all (|201lh find that the data favour a broad 
range of peak frequencies, so broad that any trend with lumi- 
nosity is blurred. They also argue for different distributions 
of break frequencies, u m , for BL Lacs and FSRQs, in the 
sense that the former objects have substantially higher val- 
ues of v m , implying that their synchrotron emission comes 
from more c ompact region s . Both points are supported by 
the study of iGiommi et al.1 (l201lh . 

The |Pe Zotti et al.1 fcooa ) model predicts 0.21 
blazars/deg 2 brighter than 35 mjy at 500 /im, i.e. ~ 28 
blazars over the H-ATLAS Phase 1 area. At the same flux 
density limit the ITucci et al.l (|201lf ) predictions are in the 
range 0.068-0.12 blazars/deg 2 , implying about 9-16 blazars 
over the Phase 1 area, the lower values referring to the 
C2Ex model, the higher to the C2Co model. 

As discussed in Section[2] in the H-ATLAS equato- 
rial fields totalling ~ 135 deg 2 we have found 9 catalogued 
blazars and 12 blazar candidates. Two of the catalogued 
blazars show signs of a substantial host galaxy contribution 
to the 500 fim flux density, so that the pure non-thermal 
emission could well be below our flux density threshold. As 
for candidate blazars, there are indications that the 500 fim 
flux density of most (8) of them is dominated or at least 
heavily contaminated by thermal dust emission (see Sec- 
tion^. 
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The contamination in the optical/near-IR translates 
into anomalous colours and may thus introduce incomplete- 
ness in the blazar selection based on diagnostics of the kind 
discussed in Section [2~Tl Moreover, the combination of ther- 
mal dust and synchrotron emission may result in very red 
sub-mm colours, so red that the H- ATLAS source extraction 
procedure, which starts from the identification of > 2.5a 
peaks at 250^m, may miss objects brighter than our 500 /im 
threshold. One such object was indeed recovered by looking 
for sign als in the 50 nm m aps at the positions of known 
blazars. iRigbv et all (|201lh estimate an incompleteness of 
12 percent for > 5a 500 /im sources. The incompleteness 
may be higher at the 4a limit adopted here and for particu- 
larly red sub-mm sources. On the other hand, any additional 
contribution to the observed flux densities biases high the 
counts compared to the case of pure synchrotron, envisaged 
by the models. 

The present analysis suggests that the second effect (ex- 
cess 500 fim flux density due to the contribution of dust in 
the host galaxy) is the dominant one. If only 4 of our blazar 
candidates have the 500 fim flux density dominated by the 
non-thermal emission, adding them to the 7 catalogued 
blazars with probably uncontaminated 500 fim flux density, 
we end up with ~ 11 blazars brighter than Ssoo/jm = 35mJy 
in our area, in good agreement with the predictions of 



iTucci et alj (|201lh C2Ex model, which is also favoured by 
other data. In any case, the 500 fim blazar counts at our 
flu x density limit are sub stantially below the predictions of 
thelDe Zotti et alj ll2005l) model. This adds to the growing 
evidence fe.g. Ipianck Collaboration et all l2011bl ) that the 
synchrotron peak frequency of most blazars occurs at lower 
frequencies than envisaged by the 'blazar sequence' scenario. 



5 CONCLUSIONS 

The sub-mm selection emphasizes blazar flavours that are 
marginal or missing in the familiar radio or X-ray selected 
samples. Two out of the 9 known blazars in the H- ATLAS 
fields (# 1 and # 6) and 8 of the 12 blazar candidates with 
S^oo^m ^ 35mJy detected by the H- ATLAS survey in its 
equatorial fields show evidence of a dust emission peak at 
sub-mm wavelengths. While we must caution that, as noted 
in Section 12.11 the H- ATLAS identification of the blazar # 1 
needs to be further checked, and that the nature of candidate 
blazars is still to be established and most of them might not 
be blazars at all, the indication that some blazar hosts are 
endowed with active star formation is supported by the fact 
that the sub-mm continuum of 5 of the 12 catalogued Planck 
ERCSC blazars detected at both 545 and 857 GHz appear 
to be dominated by thermal dust emission. This shakes the 
notion that blazar hosts are passive ellipticals. Also, some- 
what surprisingly, several of these sub-mm selected blazars 
with possibly/likely star- forming hosts are classified as flat- 
spectrum radio quasars, although quasar-type blazars are 
expected to outshine the host galaxy. 

The thermal dust contamination of the sub-mm emis- 
sion complicates the comparison of the observed blazar 
counts with model predictions. On one side the additional 
contribution to the observed flux densities leads to an over- 
estimate of the counts compared to the case of pure non- 
thermal emission, to which models refer. On the other side, 



the contamination distorts the colours so that the blazars 
may be missed by diagnostics of the kind discussed in Sec- 
tion l2.ll The present investigation however suggests that the 
surface density of blazars bri ghter than 35 mjy at 500 fim 
is lower than predicted by the De Zotti et alj (|2005h model, 
based on the blazar sequence scenario that envisages an anti- 
correlation between the synchrotron peak frequencies and 
the radio luminosity. On the other hand , such surface den - 
sity is consistent with the predictions bv lTucci et al.l (|201ll ) 
who adopted a broad distribution of synchrotron peak fre- 
quencies at all luminosities, which would have the effect of 
lowering the effective synchrotron peak frequency for the 
bright sources of interest here. 

The sub-mm selection may also bring to light blazars 
with other peculiar properties. In at least one case (blazar 
# 6, classified as a flat-spectrum radio quasar) there is ev- 
idence of the coexistence of the far-IR/sub-mm bump at- 
tributable to the host galaxy with the UV bump interpreted 
as the thermal emission from the accretion disk. Again this 
is at odds with the notion that the detection of the UV 
bump should generally be alternative to the detection of the 
emission from the host galaxy. In the case of FSRQ the lat- 
ter is outshone by the AGN emission, while in the case of 
BL Lacs the weak (or absent) line emission is indicative of 
a faint thermal emission from the disk. Thus, the accretion 
disk emission should be detectable only for FSRQs and the 
host galaxy should be visible only for BL Lacs. 

The rest-frame synchrotron peak frequencies (in terms 
of vL v ) of the 9 catalogued H-ATLAS blazars are in 
the range 12.8 ^ log f pca k,Hz ^ 14.2, with a median 
log i/poak, median, Hz — 13.4. These objects are thus mostly 
LSPs, but with an extension in the ISP region. 

At this stage all conclusions must be considered as only 
tentative because of the poor statistics. This will improve as 
soon as data for the full H-ATLAS survey, covering an area 
4 times larger, will be available. The sub-mm selected blazar 
sample will be further augmented by other Herschel surveys 
such as the Herschel Multi-tiered Ex tragalactic Survey (Her- 
MES. lHerMES Collaborationll2012[) and the Herschel Virgo 
Cluster Survey fHeVICS. lDavies et al.ll2012T l. A complemen- 
tary view of sub-mm selected blazars is being provided by 
the Planck surveys, which cover the whole sky with detec- 
tion limits more than an order of magnitude brighter than 
achieved by the H-ATLAS survey at similar wavelengths. 
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H- ATLAS IAU ID 


RA 


DEC 




>S350/jm 


S , 500fim 
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z 


BZCAT NAME 


Fermi 2FGL NAME 








roi 


r t i 
[mjy] 


[mjy] 


r t i 
[mjy] 


r t i 
[mjy] 


[arcscc] 








1 


HATLAS J083949.3+010436 


129.95576 


1.07669 


50.50 


50.68 


42.83 


443.71 


9.72 


1.123 


BZQJ0839 + 0104 


2FGL J0839.6+0059 


2 


HATLAS J090940.3+020000 


137.41804 


2.00013 


40.29 


64.49 


73.87 


305.76 


8.67 


~ 


BZB J0909 + 0200 


2FGL J0909.6+0158 


3 


HATLAS J090910.1+012135 


137.29245 


1.35986 


258.11 


337.74 


424.86 


559.64 


1.42 


1.024 


BZQ J0909 + 0121 


2FGL J0909.1+0121 


4 


HATLAS J115043.8-002355 


177.68276 


-0.39885 


37.90 


63.28 


79.22 


2803.17 


1.90 


1.976 


BZQJ1150 - 0023 




5 


HATLAS J113245. 7+003427 


173.19054 


0.57434 


65.81 


73.66 


57.25 


468.98 


1.43 


" 


SZSJ1132 + 0034 


2FGL J1132. 9+0033 


6 


HATLAS J113302.9+001545 


173.26229 


0.26256 


59.41 


48.56 


40.35 


214.96 


3.99 


1.173 


BZQJ1133 + 0015 
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HATLAS J113320.1+004054 


173.33408 


0.68185 


35.61 


43.24 


46.90 


312.49 


2.40 


1.633 


BZQJ1133 + 0040 
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HATLAS J141004.6+020306 
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114.23 


152.90 


178.05 


291.62 


0.80 




BZBJUW + 0203 




9 


HATLAS J140412. 1-001325 


211.05040 


-0.22360 


18.00 


32.00 


45.00 


516.32 


0.08 


1.217 


BZQJ1A0A - 0013 
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Table 3. BZCAT blazars found in the H- ATLAS 9h, 12 h and 15 h fields. Objects BZBJ0909+0200 and BZQJ0909+0121 are those found by Gonzalez-Nuevo et al. (2010) in the 
H- ATLAS Science Demonstration Phase field. Objects #1-8 are listed in H- ATLAS catalogues while object #9 lies near a border of the map and was excluded when making the 
catalogues. The distance between the H-ATLAS source and the FIRST low radio frequency counterpart is shown. The rcdshifts are from the BZCAT. 



ID H-ATLAS IAU ID RA DEC S 2 50 M m S350 M m Ssoo M m Sfirst Distance FIRST ID FIRST RA FIRST DEC 

[°] [°] [mjy] [mjy] [mjy] [mjy] [arcsec] [hms] [dms] 



10 HATLAS J143405. 2-002307 218.52182 -0.38549 73.40 70.02 41.95 59.92 3.01 J143405. 2-002310 14 34 05.248 -00 23 10.76 



11 HATLAS J144119. 7+002511 220.33240 0.41993 52.44 60.79 41.15 57.32 7.31 J144119.4+002506 14 41 19.412 +00 25 06.88 



12 HATLAS J140931. 9-010054 212.38320 -1.01520 43.92 38.33 35.31 63.61 7.60 J140931. 4-010053 14 09 31.465 -01 00 53.83 



13 HATLAS J144930. 2+023636 222.37620 2.61022 24.43 43.89 46.27 50.09 4.32 J144930.3+023632 14 49 30.310 +02 36 32.48 



14 HATLAS J113622. 0+004854 174.09173 0.81505 99.44 65.30 43.97 151.74 2.40 J113622.0+004851 11 36 22.041 +00 48 51.80 



15 HATLAS J085418. 9-003612 133.57879 -0.60352 56.07 68.31 48.88 25.43 7.24 J085418. 6-003619 08 54 18.690 -00 36 19.14 



u 16 HATLAS J114616. 8-001937 176.57016 -0.32712 95.44 69.55 37.28 22.17 1.11 J114616. 8-001936 11 46 16.820 -00 19 36.53 
o 

o 

8 17 HATLAS J114203. 3+005133 175.51381 0.85938 147.77 66.52 36.79 15.92 3.00 J114203.4+005136 11 42 03.440 +00 51 36.05 



> 18 HATLAS J142631. 1+014226 216.62988 1.70600 30.70 43.20 49.80 23.39 5.92 J142631.4+014231 14 26 31.400 +01 42 31.25 
w 

g 19 HATLAS J143743.4-005237 219.43098 -0.87711 53.80 49.90 47.70 35.40 8.76 J143742.8-005240 14 37 42.880 -00 52 40.64 



£ 20 HATLAS J145146. 2+010610 222.94281 1.10293 35.60 60.30 62.70 41.78 9.83 J145145.9+010619 14 51 45.950 +01 06 19.10 

w 

o 21 HATLAS J144937.4+004135 222.40611 0.69322 53.42 44.550 46.26 51.93 8.19 J144937.9+004131 14 49 37.946 +00 41 31.66 

o 

o Table 4. New blazar candidates found in the H-ATLAS 9h, 12 h and 15 h fields. Rcdshifts are not available for any of these objects, 
o 

i 
o 
o 
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3618.51 
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343.4920 


16.1481 
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18340.72 


12657.00 


0.859 


BZQJ2253+1608 




187.2750 


2.0458 


6937.50 
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4065.28 


54991.00 
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BZQJ1229+0203 
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-5.7894 


4701.58 
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9711.00 
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-36.0869 
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1.417 


BZQJ0403-3605 




84.7120 


-44.0853 


4065.40 


3826.54 


3057.01 


3729.00 


0.892 


BZBJ0538-4405 




140.3960 


62.2628 


1234.73 


1495.28 


3041.21 


946.00 


1.446 


BZQJ0921+6215 




130.1979 


13.2066 


2801.51 


3263.22 


5745.42 


2614.00 


0.681 


BZUJ0840+1312 


2FGLJ0840.7+1310 



Table 5. BZCAT blazars with a P/ancfc-ERCSC counterpart at both 545 and 857 GHz (see text). The 600 GHz flux densities were 
computed interpolating between measurements at 545 and 857 GHz. The rcdshifts arc from the BZCAT. 
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